String inspired explanation for the super-acceleration of our universe 
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We investigate the effect of the bullc content in the general Gauss-Bonnet braneworld on the 
evolution of the universe. We find that the Gauss-Bonnet term and the combination of the dark 
radiation and the matter content of the bulk play a crucial role in the universe evolution. We show 
that our model can describe the super-acceleration of our universe with the equation of state of the 
effective dark energy in agreement with observations. 
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I ■ I. INTRODUCTION 

^ ■ 

' One of the most dramatic discoveries of the modern cosmology in the past decade is that our universe is currently 
, accelerating fl'| . A component that causes the accelerated expansion of the universe is referred to as "dark energy" . 
Although the nature of such dark energy is still speculative, an overwhelming flood of papers has appeared which 
r~~i attempt to describe it by devising a great variety of models. Among them are cosmological constant, exotic fields such 
'"^^ as phantom or quintessence, modified gravity, etc, see y for a recent review. Available models of dark energy differ 
' ^ in the value and variation of the equation of state parameter w during the evolution of the universe. The cosmological 
constant with w = located at a central position amongdark energy models both in theoretical investigation 

r-| ] and in data analysis [1| . In quintessence Q , Chaplygin gas [1] and holographic dark energy models [3] , w always 
stays bigger than —1. The phantom models of dark energy have w < — 1 [8|. However, following the more accurate 
data analysis, a more dramatic result appears showing that the time varying dark energy gives a better fit than a 
cosmological constant and in particular, w can cross —1 around z = 0.2 from above to below [9|. Although the galaxy 
cluster gas mass fraction data do not support the time- varying w (lo| , theoretical attempts towards the understanding 
\^ ' of the w crossing —1 phenomenon have been started. Some dark energy models, such as the one containing a negative 
\^ • kinetic scalar field and a normal scalar field fll|, or a single scalar field model and interacting holographic dark 
' energy models [l3| have been constructed to gain insight into the occurrence of the transition of the dark energy 
. equation of state and the mechanism behind this transition. Other studies on the w = — I crossing have been carried 
O ■ out in [III. 

C — , Independent of the challenge we deal with the dark energy puzzle, in recent years, theories of large extra dimensions, 
' in which the observed universe is realized as a brane embedded in a higher dimensional spacetime, have received a lot 
^ ' of interest. According to the braneworld scenario the standard model of particle fields are confined on the brane while, 
in contrast, the gravity is free to propagate in the whole spacetime. In these theories the cosmological evolution on the 
brane is described by an effective Friedmann equation that incorporates non-trivially with the effects of the bulk onto 
the brane. An interesting consequence of the braneworld scenario is that it allows the presence of five-dimensional 
matter which can propagate in the bulk space and may interact with the matter content in the braneworld. It has 
been shown that such an interaction can alter the profile of the cosmic expansion and lead to a behavior that would 
resemble the dark energy. The cosmic evolution of the braneworld with energy exchange between brane and bulk has 
been studied in different setups [iB] [I3| O • these models, due to the energy exchange between the bulk and 
the brane, the usual energy conservation law on the brane is broken and consequently it was found that the equation 
of state of the effective dark energy may experience the transition behavior (see e.g pTi. [isj). 

On the other hand, in string theory, in addition to the Einstein action, some higher derivative curvature terms have 
been included to derive the gravity. In order to obtain a ghost-free theory, the combination of quadratic terms called 
Gauss-Bonnet term is usually employed as curvature corrections to the Einstein-Hilbert action [19]. From a geometric 
point of view, the combination of the Einstein-Hilbert and Gauss-Bonnet term constitutes, for 5D spacctimcs, the 
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most general Lagrangian to produce second-order field equations [20|. The Gauss-Bonnet correction significantly 
changes the bulk field equations and leads to modifications in the braneworld Friedmann equations. Therefore, the 
study of the effects of the Gauss-Bonnet correction term on the evolution of the universe in the braneworld scenario 
is well motivated. Influences of the Gauss-Bonnet correction on the DGP braneworld have been studied in [21I, . 

The purpose of the present work is to investigate the effects of the bulk content in the general Gauss-Bonnet 
braneworld on the evolution of the universe. Although the effects of the Gauss-Bonnet correction term on the late 
time universe is small, we will see that it still plays an important role in the cosmic evolution. Besides we will 
show that the combination of the dark radiation term and the matter content of the bulk plays the role of the dark 
ener gy o n the brane and influences the evolution of the universe. In our model, in contrast to the previous models 
([isl llGl [TtI . nil), we do not need to break down the standard energy momentum conservation law on the brane, 
although our model can allow such assumption if one is interested. We will show that by suitably choosing model 
parameters, our model can exhibit accelerated expansion of the universe. In addition, we will present a profile of the 
w crossing —1 phenomenon which is in good agreement with observations. 

The paper is organized as follows. In SectionHIl we present a braneworld model to describe the accelerated expansion 
and the effective equation of state of dark energy in the presence of the Gauss-Bonnet correction term in the bulk. In 
Section [ml we study the cosmological consequences of the model and in particular, its effect on the evolution of the 
universe. The last section is devoted to conclusions and discussions. 



II. THE MODEL 



The theory we are considering is five-dimensional and has an action of the form 

5 = ^ y d'x^ (i? - 2A + aCGB) + J d^x^CZik + j d^x^g{CZane - <y). (1) 
where A < is the bulk cosmological constant and Cg b is the Gauss-Bonnet correction term 

-Cgb = -R^ — AR^^ Rab + R^^'^^ Rabcd ■ (2) 

Here g and g are the bulk and brane metrics, respectively. R, Rab, and Rabcd are the scalar curvature and Ricci 
and Riemann tensors, respectively. Throughout this paper we choose the unit so that = 1 as the gravitational 
constant in five dimension. We have also included arbitrary matter content both in the bulk and on the brane through 
^buik ^^'^ ^brane respectively, and a is the positive brane tension. The field equations can be obtained by varying the 
action ([T]) with respect to the bulk metric gAB- The result is 

Gab + J^gAB + '2aHAB = Tab, (3) 
where Hab is the second-order Lovelock tensor 

Hab — RRab — '2'Ra'^ Rbc ~ '2,R'~^^ Racbd 
+Ra^^^ Rbcde — j9ab^gb ■ 

For convenience and without loss of generality, we can choose the extra-dimensional coordinate y such that the brane 
is located at y = and bulk has Z2 symmetry. We are interested in the cosmological solution with a metric 

^ -n^(t,y)dt^ +a'^{t,y)-f^jdx'dx^ + b'^{t,y)dy'^, (4) 

where 7^ is a maximally symmetric 3-dimensional metric for the surface (t=const., ?/=const.), whose spatial curvature 
is parameterized by k = -1, 0, 1. The metric coefficients n and b are chosen so that, n{t, 0) = 1 and b{t, 0) — 1, where t 
is cosmic time on the brane. The total energy-momentum tensor has bulk and brane components and can be written 
as 

Tab ~ Tab Wane +Tab U +Tab \buik ■ (5) 

The first and the second terms are the contribution from the energy-momentum tensor of the matter field confined 
to the brane and the brane tension 

Tji \brane = diag(-p, p, p, p, 0) , (6) 

Til, = diag(-a,-a,-a,-a,0)^, (7) 
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where p, and p, being the energy density and pressure on the brane, respectively. In addition we assume an energy- 
momentum tensor for the bulk content of the form 



rpA 



B \bulk - 



J- f 



n'' rpO 
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The quantities which are of interest here are T\ and T\, as these two enter the cosmological equations of motion. In 
fact, T\ is the term responsible for energy exchange between the brane and the bulk. Integrating the (00) component 
of the field equations ^ across the brane and imposing Z2 symmetry, we have the jump across the brane [23| 



l + Aa\ H 
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where 2a'^ ~ ~2a'_ is the discontinuity of the first derivative. H — ciq/ qq is the Hubble parameter on the brane. Eq. 
([9]) is a cubic equation for the discontinuity a^/ao, which has only one real solution, the other two being complex. 
Therefore, if we require our cosmological equations to have the right a — > limit we are left with only one solution. 
However, this real root is too lengthy and complicated to present here. Since we are interested to study the effect of 
the Gauss-Bonnet correction term on the evolution of the universe in the late time so it is reasonable to choose the 
Gauss-Bonnet coupling constant a to be small, namely < a < 1. Using this fact we can expand the real solution 
for a'j^/a^ versus a powers. The result for fc = up to order a is 
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In a similar way, integrating the (ij) component of the field equations Q across the brane and imposing Z2 symmetry, 
we can obtain the discontinuity in the metric function rt^ /no, which for k = can be written up to 0(q;) in the following 
form 



n'. 



no 



(2p + 3p-a) + -(-2i?2 (2p4-3p- cr) 
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where dots denote time derivatives and primes denote derivatives with respect to y. At this point we find it convenient 



to absorb the brane tension cr in p and p with the replacement p 
conditions pI7|) and ([TT|) can be simplified 



(7 ^ p and p — a ^ p. Therefore the junction 
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Substituting the junction conditions and (|13p into the (55) and (05) components of the field equations we 
obtain the modified Friedmann equation and the semi-conservation law on the brane (up to order a) 



1-^P(2P 
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1 + 4a Lff + — 
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and 



p + ZH{p + p) = -T, T=2TH 
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We shall assume an equation of state p 
Therefore we have 



wp to hold between the energy density and pressure of matter on the brane. 
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One can easily check that in the hmit a ^ 0, Eqs . ((T0l) - ((T7)) reduce to the corresponding equations of the braneworld 
model without Gauss-Bonnet correction term [l^ . 

Remarkably, we can show that the Friedmann equation (fT6|) is equivalent to the following equations 



with X satisfying 
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Using the definition for T in Eq. (jlSp . the latter equation up to order a can be written as 
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(18) 



(19) 



(20) 



Eq. (jl8[) is the modified Friedmann equation describing cosmological evolution on the branc. The auxiliary field 
X incorporates non-trivial contributions of dark energy which differ from the standard matter fields confined to the 
brane. The bulk matter contributes to the energy content of the brane through the bulk pressure terms T\ that 
appear in the right hand side of the Friedmann equation. In addition, the bulk matter contributes to the energy 
conservation equation (|15|) through T\ which is responsible for the energy exchange between the brane and bulk. 
The functions T\ and T\ are functions of time corresponding to their values on the brane. The energy-momentum 
conservation \7 = cannot fully determine T\ and T\ and a particular model of the bulk matter is required 
[ri |. In the limit a ^ 0, Eqs. p8)) and ([20]) reduce to ( after replacement p ^ p + a) 
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If we invoke the usual definitation (3 =\/ 36, A 



7^5 
6 

I (A + ctV6)/6 and 7 = cr/3, we get (k^ = l) 



(21) 
(22) 
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which is noting, but the general set of the equations in RS II braneworld model with bulk matter content plus 
brane-bulk energy exchange (see for example [18]). 

Returning to the general Friedman equation (|18p with Gauss-Bonnet correction term, we can show that this equation 
has the solution for H of the form 



H' 



1 

Aa 



36 108a 



[729 -I- 12ap2(27 + ap^) + 972a(6x + A - 2T|)] 



1/2 



The upper solution (+) has correct a limit. Indeed, if we expand this solution verses a we get (up to 0(a)) 
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and Eqs. ((151) ^^^^ pOl) become 

p + ?.Hp{l+uj) = -2T% 
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(24) 

) 

(25) 

(26) 
(27) 



Therefore, until now we have obtained the set of equations describing the dynamics of our universe (Eqs. (|25 p -(|27 p ) 
in the general Gauss-Bonnet braneworld with both bulk matter content and bulk-brane energy exchange provided 
that the Gauss-Bonnet coupling constant a is chosen sufficiently small. It is worth noting that although a is small, 
it has a dramatic effect on the dynamic behavior of the cosmic evolution. Besides the appearance of the term on 
the right hand side of Eq. (I25p shows that in high energy scale the Gauss-Bonnet correction term plays an important 
role. 
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Figure 2: Evolution of Weff(^) versus z for < v < 0.17. a = 0.01 (bold line), a = 0.1 (continuous line), and 
line). 
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Figure 3: Evolution of q{z) versus z ior < v < 0.17. a = 0.01 (bold line), a = 0.1 (continuous line), and a = 0.9 (dashed 
line) . 
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Figure 4: Evolution of WeH{z) versus z for 0.18 < v < 0.34. a = 0.01 (bold line), a = 0.1 (continuous line), and a = 0.9 
(dashed lined). 
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Figure 5: Evolution of q{z) versus z for 0.18 < v < 0.34. a = 0.01 (bold line), a = 0.1 (continuous line), and a = 0.9 (dashed 
line) . 
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Figure 6: The parameter space of the function a{i',A) for < z/ < 0.17. ^ = 18 (bold line), j4 = 30 (continuous line), and 
A = 50 (dashed line). 
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Figure 8: The parameter space of the function a{i',A) for 0.18 < v < 0.34. A = —6 (bold Une), A = —10 (continuous line), 
and A = —20 (dashed Une). 
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Figure 9: The parameter space of the function C{v,A) for 0.18 < u < 0.34. A = —6 (bold line), A — —10 (continuous line), 
and A = —20 (dashed line). 



III. COSMOLOGICAL CONSEQUENCES 



In this section we are going to explore some cosmological consequences of our model. To do this, first we separate 
back the matter energy density and the brane tension as usual form with the replacement p p + a. Therefore Eqs. 
13 and (HZI) become 



H' = 
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We are interested in the scenarios where the energy density of the brane is much lower than the brane tension, namely 
/9 <C cr. Assuming the Randall-Sundrum fine-tuning A + = holds on the brane and defining the parameter 
7 = (t/36, Eqs. and can be simplified in the following form 



9 T\ a 
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Now, one may adopt several strategies to find solutions of Eqs. (pS)) . ([50)1 and (PT|) . For example, one may take 
a suitable ansatz for the time dependent functions T\ and T\ and using Eq. (|3ip to find the function x- Then 
substitute x, r']^ and T^g into Eq. ^ one can try to obtain p, and finally one may find Hubble parameter H through 



10 



Eq. (poll . In the following we are interested in the case in which the energy momentum conservation law on the 
brane holds, which is usually assumed in the braneworld scenarios. Indeed, we want to consider the effect of the 
bulk content on the evolution of the universe without brane-bulk energy exchange, therefore we set — 0. The 
case with brane-bulk energy exchange in the general Gauss-Bonnet braneworld will be addressed elsewhere. It was 
argued that the energy exchange between the bulk and brane T\ will lead to the effective dark energy equation of 
state crossing —1 [16,17]. Here we will show that without the energy exchange, the effect of and the combined Tg 
and the Gauss-Bonnet correction have the same role. 

Inserting the condition T\ = in Eq. it reduces to p^ZH piX^uj) — 0. This equation has well known solution 



P = Paa 



'3(1-|-Ml) 



, where po is the present matter density of the universe and we have omitted the "o" subscript from 



the scale factor on the brane for simplicity. Then, consider a general ansatz T\ = Day for the bulk pressure [15 
where D and v are two arbitrary constants, one can easily check that Eq. (pij) has a solution of the form 



X^Ca 



(32) 



where C is a constant usually referred to as dark radiation term and B = D{i' A- 2)/{iv -t- 12). Finally, inserting p 
and X ij^to Eq. pop . we can rewrite it in the standard form 
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37/47r is the 4-dimensional Newtonian constant and p^s represents the effective dark energy density on 
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(34) 



where A = —2D/ {Si' + 12). The equation of state parameter of the effective dark energy on the brane can be defined 
by m 



Weft 
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1 dlnSH^ 
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(35) 



where SH^ — {H'^ / Hq) — flmO. ^ accounts for terms in the Friedmann equation except the brane matter with equation 
of state Wm = 0. Now, if we use the redshift parameter 1 + z = as our variable, we can easily show that 
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The corresponding late time deceleration parameter can be written 

(1 - (1 + 3Weff(2))] , 



(36) 



(37) 



where il,n = ^mo (1 + ^) is all part of the matter on the braneworld and we take its present value as f2mO = 0.28±0.02. 
In the rest of the paper, we will obtain constraints on the parameters such as C, A, v, a and a in our model. Indeed, 
we want to show that under what parameter space constraints our model can describe the accelerated expansion of 
the universe with the equation of state of the effective dark energy Wcff crossing —1, as suggested by observations. 



A. Special case with a = 



Let us begin with the special case, in which the Gauss-Bonnet coupling constant a is equal to zero. In this case we 
have the usual Randall- Sundrum II braneworld model and Eq. ((36l) reduces to 



1 Uc{i+zr^^_AA 

■Wctfiz) = — 1 H r-j . (38) 

^ ' 3 I C{l + zy+^ + A i ^ ' 
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Therefore, we are left with three parameters C, A, v, and two of them are independent. Requiring that at the present 
moment Weff(z = 0) = —1.06 and w crossed —1 around z — 0.2 as indicated by extensive analysis of observational 
data [§], we can obtain 

C = 0.039A, 1^ = 0.34, A^A. (39) 

For these value of parameters and VlmQ = 0.28, from Eq. ([37]) we have q{z = 0] = —0.64 and in addition q{z) crosses 
around z = 0.33 which is in good agreement with recent observational data flKl^]. In figure [T] we plot Wcff{z) and 
q{z) for the above value of the parameters versus redshift parameter z. 

B. General case with a 

Next, we consider the general Gauss-Bonnet braneworld with bulk matter content. In this case we have five 
parameters only four of which are independent. Considering that the value of a does not affect the general profile 
of our model and further according to the Randall-Sundrum fine-tuning relation it should be small, we first fix 
a — 10~^. Thus, we have now four parameters and among them three are independent. Numerical calculations 
show that the functions WeS and q are well behaved for z > 0, provided that < i' < 0.34. Employing the present 
value of the equation of state parameter of dark energy Wcsiz = 0) = —1.06 and the moment it crossed —1, namely 
^^ctfiz — 0.2) — —1, we get 

a = a{A,iy), C ^ 0.12 Ai^ {1.2)- " , A = A. (40) 

If we impose the condition < a < 1 which was used in deriving our equations, we can get constraint on the 
free parameter A. In numerical calculations we find that for < < 0.17 we should have A > 17.97, while for 
0.18 < ly < 0.34 we should have A < —5.23 to satisfy the condition on a. In figures [5] and [3] we plot Wcb{z) and q{z) 
for < ^ < 0.17 versus redshift parameter z for different value of the Gauss-Bonnet coupling constant a. From these 
figures we observe that at large z, the Wcb{z) increases with the increase of a, while q(z) decreases with the increase 
of a. This qualitative behavior is quite opposite when 0.18 < < 0.34 as one can see from figures |4] and El Finally 
we plot in figures the parameter space for the functions a — a{A, v) and C = C{A^ v). We find that in the case 
< < 0.17, a and C increase with the increase of v while, in contrast, for 0.18 < v < 0.34, a and C decrease with 
the increase of v. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this work we have generalized the Randall-Sundrum II braneworld with both bulk matter content and bulk-brane 
energy exchange by adding the Gauss-Bonnet curvature correction term in the bulk action. We have investigated the 
effects of the bulk content in the general Gauss-Bonnet braneworld on the evolution of the universe and found that 
although the effect of the Gauss-Bonnet correction term in the late time universe is small, it still plays an important 
role in the universe evolution. 

In contrast to the previous models ([l^, [l^ 113, IHDj in our study we kept the energy momentum conservation 
law on the brane as usual and found that the combination of the dark radiation term and the matter content of 
the bulk can play the role of the dark energy on the brane and influence the evolution of the universe. By suitably 
choosing parameter space in our model, we can describe the super-acceleration of our universe with the behavior of 
the effective dark energy equation of state in agreement with observations. In [2^, it was argued in a Gauss-Bonnet 
brane world with induced gravity that the Gauss-Bonnet term and the mass parameter in the bulk play a crucial 
role in the evolution of the universe. Here in our general model, we confirmed their argument. It is easy to see 
from Eqs.(25)-(27) that the Gauss-Bonnet correction influences the dynamics of our universe, especially in the early 
universe with high energy scale. Phenomenon on the Gauss-Bonnet role has been disclosed in Figs. 2-5. We observed 
that although the Gauss-Bonnet effect is not clear at the present moment, it influences the universe evolution in the 
past and was more important in the earlier period. 

In this work we just restricted our numerical fitting to limited observational data. Giving the wide range of 
cosmological data available, in the future we expect to further constrain our model parameter space and test the 
viability of our model. 
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